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ABSTRACT 

This paper reviews recent progress of lattice gauge theory determinations of 
semileptonic B and D decay form factors. These determinations are important 
in extracting the remaining CKM matrix elements. 



1 Introduction 

A large scale experimental and theoretical effort is underway to overconstrain 
the CKM matrix and uncover any internal inconsistencies revealing new physics. 
Uncertainties in determinations of the CKM elements arise not only from ex- 
periment but also, in no small part, from the theoretical calculations needed to 
account for hadronic QCD effects. Within the next few years the theoretical 
uncertainties must be reduced to the few percent level in order not to dominate 
uncertainties from experiment. 



B decay and mixing processes are most suitable for extracting \ Vub\, \Vcb\, 
\Vtd\, and \Vts\', a large fraction of the CKM matrix. The QCD coupling at the 
relevant hadronic scales is large and thus a non perturbative method is needed. 
Lattice QCD holds the most promise in being able to provide the required 
hadronic factors with the required accuracy. Recent lattice QCD results give 
confidence that the required precision calculations are now possible. 

This paper reviews the current status of lattice QCD calculations for 
the semileptonic decays of the B mesons needed to determine Vub and Vcb- 
Additionally, semileptonic D meson decay calculations, suitable as rigorous 
checks of lattice methods, are discussed. 

2 Lattice QCD 

Lattice QCD (see, e.g. ^1) involves the use of a mathematical 'trick' where 
spacetime is discretized into a finite lattice. Quarks live on the lattice points 
and gluons live on the links between the points. This formalism regularizes 
QCD by providing a momentum cut-off: no momenta greater than n/a can 
propagate where a, typically ~ O.lfm, is the lattice spacing. The Feynman 
path integral becomes an ordinary integral over a finite number of degrees 
of freedom, and can be computed numerically on a computer. Continuum 
QCD results can be obtained by taking the lattice spacing to zero, provided 
that the matching factors, or differences between the continuum and lattice 
renormalization schemes, have been taken into account. In heavy-light physics 
these matching factors are usually determined by comparing a perturbative 
continuum calculation with the corresponding perturbative calculation on the 
lattice. 

Unfortunately, the numerical integrations corresponding to the Feynman 
path integral are extremely computationally expensive. Even with the use of 
efficient Monte Carlo methods approximations must be done in order to obtain 
results with the computational technology of today. A dramatic saving can be 
made by ignoring closed quark loops in the vacuum, and the vast majority of 
lattice calculations have been done in this so called quenched approximation. 
The use of this incorrect theory, however, leads to systematical errors at the 
10-20% level. Unquenched calculations must be done in order to achieve the 
above precision results. 

In unquenched calculations, when vacuum or dynamical quarks are in- 



eluded, the expense of the simulation increases dramatically with decreasing 
dynamical quark mass, meaning that in practice the light dynamical quarks are 
included with masses greater than their physical masses. If light enough, how- 
ever, extrapolations can be done to the correct physical masses with the help 
of 'chiral perturbation theory' (see, e.g. EJ): an effective theory involving ex- 
pansions around the massless limit. Up until recently, unquenched simulations 
have not been able to reach this 'chiral regime'. 

Now for the first time, simulations have been done with dynamical quarks 
light enough to allow the agreement at the 3% level of theory with experiment 
for a variety of (simply calculable) quantities IS*. These simulations have been 
possible due to the combination of ever increasing computing power and the 
emergence of a better understanding about the properties of quarks on the 
lattice, which has lead to the use of the so called improved sto^^ered formulation. 
This was used by the MILC collaboration (seeS ^ and references therein) to 
create ensembles of 'configurations' (snapshots of the QCD vacuum on the 
lattice) which are then used to 'measure' required physical quantities such as 
those above. 

3 Heavy Quarks on the Lattice 

Naive discretization of heavy quark fields leads to large 0{amQ) discretization 
errors due to large heavy quark mass mq. However, one can use effective theo- 
ries which take advantage of the fact that the heavy quarks typically have low 
velocities within the hadron, and are therefore somewhat non-relativistic. This 
also often results in simplifications reducing simulation time and allowing high 
statistics. The calculations reviewed in this paper incorporate two alternative 
heavy quark methods. 

Non Relativistic QCD (NRQCD) involves an expansion of the QCD 
Lagrangian in powers of l/mq. This is very useful for h quarks but not so 
appropriate for c quarks. 

The Fermilab methodic (although more complicated than NRQCD) is 
very appropriate for c quarks as it incorporates smooth transitions from rela- 
tivistic light quarks to non relativistic heavy quarks. 



4 Semileptonic Decays 



Recently the first fully unquenched results for the semileptonic decay form 
factors have appeared and, although preliminary, are very promising. 

The matrix element for the decay of a heavy B ov D meson to a pion is 
given by 



{tt{p^)\V^\H{ph)) = f+{q^) 



■,2 



q 

= V2^ {f\\iE^)v^ + f^iEr,)p''J (1) 

where u*" — p^/Mn and p^ = pf^ — E-^v^^ become (1, 0) and (0,p^) respectively 
in the heavy meson rest frame. 

The alternative and /_L(£'7r) form factors have been introduced 

because they are more appropriate for lattice calculations and the associated 
chiral perturbation theory formulae are usually given in terms of E-^ , the pion 
energy in the heavy meson rest frame. It is straightforward to interchange 
between these two conventions. 

Unfortunately, experimental results are limited to the small region 
whereas, for B mesons, lattice calculations are most reliable for small recoil 
(large q^). This is because lattice calculations currently work in the B meson 
rest frame and large recoil would give the pion large momenta introducing large 
0{a?p'^) discretization errors and large statistical errors. However, the lattice 
community is excited about the development of Moving NRQCD Bl which 
allows the momentum to be shared between the B and the tt. This will allow 
calculations over the full q^ range providing excellent overlap with experiment 
and should be ready to implement within the next year. 

In the mean time it is necessary to use a model in order to extrapolate the 
lattice results to the low q^ region. It will, however, be seen that results using 
the Becirevic Kaidalov (BK) parameterization ISl are encouraging. The BK 
parameterization satisfies Heavy Quark Effective Theory scaling laws, the fact 
that /+ must have a pole at q^ = Mg, and the necessary condition /+(0) = 
/o(0). 



lin[solid], xPT[dashed], stag_xPT[dotted] fits (4 points) 




0.01 0.02 0.03 



Figure 1: Chiral extrapolations for ff^'" with Fcrmilab heavy quarks l**^. 
Sohd: linear; dashed: chiral perturbation theory without staggered effects; 
dotted: full staggered chiral perturbation theory. 



5 Semileptonic D Decay Results 

The CKM elements \Vcs\ and \Vcd\ are known more accurately than \Vub\ and 
the CLEO-c program aims to further improve this accuracy Semileptonic 
D decays are thus very suitable processes for testing lattice calculations. 

Using the MILC ensembles and Fermilab heavy quarks, Okamoto et al 1^ 
have calculated the D ^ tt (and similarly D — > K) form factors. Lattice 
determinations of the matrix element Eq. ^ were done for several and 
light quark mass m^. For each ruq, the BK parameterization was used to 
interpolate to common E^r values. Then, for each E^^, chiral extrapolations were 
done to obtain results at the physical light quark masses. Chiral perturbation 
theory is used to give the appropriate extrapolation function, dependent on 
the lattice action used. Figure ^ shows chiral extrapolations for f±. The 
correct (full staggered) chiral extrapolation is shown along with linear and non- 
staggered functions, pv is a matching factor between the lattice and continuum 
renormalization schemes which is 1 at tree-level. 

Figure 121 shows results for linearly^ chirally extrapolated form factors as 

-'^The analysis for the full staggered chirally extrapolated form factors is 
underway. 




a function of for both D ^ tt and D ^ K. In both cases good agreement 
with experiment at = is seen, although the CLEO-c results will provide a 
more stringent test over the full range. 

There exists one set of MILC ensembles with a relatively large lattice 
spacing (the 'coarse' set) and there similarly exists a 'fine' set. This work has 
so far only been done on the coarse ensembles and must be repeated on the 
fine ensembles to check for lattice spacing dependence. 



6 B ^ nlv Results 



The MILC coarse ensembles were again used for B nlv calculations. Both 
NRQCD and Fermilab heavy quarks have been used. 

In similar fashion to the D decay analysis, the form factor results were 
interpolated to common i?^ values and chiral extrapolations were performed, 
but only linearly so far. Figure O shows these extrapolations for /|| and f± 
for the NRQCD case 1^ . The full staggered chiral function has recently been 
determined and must now be incorporated into this analysis. The BK param- 
eterization was then used to extrapolate to the low region, as can be seen 




Figure 3: Linear chiral extrapolations for /|| and ff^^ with NRQCD heavy 
quarks 1^ . 



on the left hand side of figure 01 The data fits the model excellently with /+ 
exhibiting the expected pole at — M^, and with /o consistent with the soft 
pion relation /o(Af^) = ]bI I-n- Although /q is not needed for the decay rate, 
it has relatively small statistical errors and its inclusion in the fit is very useful 
in constraining /+. These results include one- loop matching but OiXj arrih) 
currents have still to be included, and again this work must be repeated on the 
fine ensembles. This plot includes old quenched results for comparison, some 
of which have had their errors removed for clarity. 

An equivalent plot of results with Fermilab quarks is shown on the 
right hand side of figure 01 again comparing with old quenched results. In this 
case again only the coarse ensembles have been used, only tree level matching 
has been done, and rnj, is not tuned well. These issues are being addressed. 

From their NRQCD form factor results, Shigemitsu et at have esti- 
mated a result for \Vub\ by integrating 

— —-^pl\U{q')\' (2) 



where the CLEO branching fraction was used to get F. The preliminary 




Figure 4: BK parameterization fits to f"^'" and f^^^ with NRQCD heavy 
quarks^J* (left), where the burst shows the soft pion result, and with Fermilab 
heavy quarks (right). Both plots include old quenched results. 



results from both the full and high q ranges are 

r 3.86(35) (62) X 10-3 Q < q^ < ql,, 

\Vub\ = <^ (3) 
[ 3.52(70)(42) X lO^^ IQGeV"^ < q^ 

where the errors are experimental and lattice respectively, and are both ten- 
tative. It is encouraging that these results agree with each other, and are 
consistent with the current inclusive B decay determinations . 

It is hoped that the lattice errors for this quantity will be at around the 
10-13% level when this analysis is complete. The main sources of error are 
uncertainties in chiral extrapolations, continuum extrapolations and matching. 
Estimates have been made as to how the magnitude of the overall error will 
reduce with future calculations . The next generation of machines, being 
built just now, should allow simulations where or mi is halved, shrinking the 
error to the 5.5-6.5% level, assuming that 2-loop matching has been performed. 
Looking further ahead, if both and mi could be halved it is hoped that 4-5% 
precision will be possible, again assuming 2-loop matching. The timescale for 
this, however, is not known. 



7 B ^ D*lv Results 



The differential B D*lv decay width is given by 

fir 

— cx iKbl^s^ij.H (4) 

where tu = v' .v and v and v' are the B and D four- velocities respectively. In 
order to extract \Vcb\, the form factor at zero recoil J-^l) must be determined. 
Since B and D mesons can both be considered heavy, heavy quark symmetry 
can be exploited. The errors then scale with 1 — ^{1) instead of because 
in the infinitely heavy quark limit J-{1) — 1 . 

The best lattice determination thus far, which is in the quenched approx- 
imation and uses Fermilab heavy quarks, is given 

where the errors are from statistics, matching, lattice spacing dependence, chi- 
ral extrapolation, and quenching respectively. 

With reference to the previous section, it is hoped that this total error 
will be reduced from 4% to around the 2% level with the next generation of 
machines and then to as low as 1% with the next again generation. At this 
level, it will be important to compute the slope and curvature of ^{uj), and 
Moving NRQCD will again be of great help in achieving this. 



8 Conclusions 



Unquenched lattice gauge theory calculations are appearing, and have already 
made an impact. The lattice community is confident that such calculations 
can now be done to obtain the quantities important for extracting the CKM 
elements, including the semileptonic B decay form factors. Comparison of D 
semileptonic decay lattice results with precise CLEO-c data should enhance 
this confidence. 

The calculations are still at a preliminary stage, but good understanding 
exists on, and plans are in place to address, all sources of error. Precision 
results are likely to appear within the next few years. 
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